These data are made available as part of paper: S. J. Gonzalez-Roji, J. Saenz, G. Ibarra-Berastegi, J. Diaz de Argandona (2018) Moisture balance over the Iberian Peninsula according to a regional climate model, Journal of Geophysical Research 123(2):708-729. The dataset holds selected postprocessed files that allow to reproduce all the results in the paper. An analysis of the atmospheric branch of the hydrological cycle by means of a 15 km resolution numerical integration performed using Weather Research and Forecasting (WRF) nested in ERA Interim is presented. Two WRF experiments covering the period 2010–2014 were prepared. The first one (N) was configured as in standard numerical downscaling experiments. The second one (D), with the same parameterizations, included a step of 3DVAR data assimilation every 6 h. Experiment D: -windowed-D-AI: Analysis increments (T2, Q2), 2010-2014 for the D experiment -windowed-wavacurip4paper-D: Netcdf files corresponding to experiment D (using 3DVAR assimilation) and years 2010-2014 holding 3-hourly outputs from the run after selecting the window covering the Iberian Peninsula. Experiment N: -windowed-wavacurip4paper-N: Netcdf files corresponding to experiment N (no data assimilation) and years 2010-2014 holding 3-hourly outputs from the run after selecting the window covering the Iberian Peninsula

OceanRAIN version 1.0, OceanRAIN-M - Number Concentration Particle Size Distribution and Precipitation Microphysics, 37 along-track parameters plus 128 size bins for 8 ships, 692.000 precipitation minutes in total, temporally discontinuous data for each ship, 1-minute-resolution

OceanRAIN version 1.0, OceanRAIN-W - Water cycle components, 73 along-track parameters for 8 ships, 6.83 million minutes in total, temporally continuous data for each ship, 1-minute-resolution

OceanRAIN version 1.0, OceanRAIN-R - ODM470 Raw number count Particle Size Distribution and Precipitation Microphysics, 37 along-track parameters plus 128 size bins for 8 ships, 692.000 precipitation minutes in total, temporally discontinuous data for each ship, 1-minute-resolution

These data are made available as part of paper: S. J. Gonzalez-Roji, J. Saenz, G. Ibarra-Berastegi, J. Diaz de Argandona (2018) Moisture balance over the Iberian Peninsula according to a regional climate model, Journal of Geophysical Research 123(2):708-729. The dataset holds selected postprocessed files that allow to reproduce all the results in the paper. An analysis of the atmospheric branch of the hydrological cycle by means of a 15 km resolution numerical integration performed using Weather Research and Forecasting (WRF) nested in ERA Interim is presented. Two WRF experiments covering the period 2010–2014 were prepared. The first one (N) was configured as in standard numerical downscaling experiments. The second one (D), with the same parameterizations, included a step of 3DVAR data assimilation every 6 h. Experiment D: -windowed-D-AI: Analysis increments (T2, Q2), 2010-2014 for the D experiment -windowed-wavacurip4paper-D: Netcdf files corresponding to experiment D (using 3DVAR assimilation) and years 2010-2014 holding 3-hourly outputs from the run after selecting the window covering the Iberian Peninsula. Experiment N: -windowed-wavacurip4paper-N: Netcdf files corresponding to experiment N (no data assimilation) and years 2010-2014 holding 3-hourly outputs from the run after selecting the window covering the Iberian Peninsula

The climatological dataset was produced using the Weather and Research Forecasting (WRF) model, version 4.2.2, configured with two nested domains at 10 km (D1) and 2 km (D2) horizontal grid spacing. It covers most of the South Island of New Zealand and is centered over Brewster Glacier in the Southern Alps. The model was forced every three hours by ERA5 reanalysis data at its outer lateral boundaries. The dataset spans the period of 1 January 2005 to 31 December 2020, providing daily output in the outer domain (D1) and 3-hourly output in the innermost domain (D2). The data provided here are a selection of daily averages from the inner WRF domain (D2; 2-km grid spacing). They are distributed among three different file types containing 4-dimensional, 3-dimensional and time-invariant output variables, respectively. For the 4-dimensional fields, perturbation and base-state atmospheric pressure (WRF variables P and PB) and geopotential (PH and PHB) were combined to produce full model fields (PRES and GEOPT). Perturbation potential temperature (T) was converted to total potential temperature (THETA). Wind vectors (U,V, and W) were converted to mass points and rotated to earth coordinates. ------- Acknowledgements: The modeling and related research was supported by the German Research Foundation (DFG) grant no. 453305163. The authors gratefully acknowledge the scientific support and HPC resources provided by the Erlangen National High Performance Computing Center (NHR@FAU) of the Friedrich-Alexander-Universität Erlangen-Nürnberg (FAU) under the NHR project b128dc / ATMOS ("Numerical atmospheric modeling for the attribution of climate change and for model improvement"). NHR funding is provided by federal and Bavarian state authorities. NHR@FAU hardware is partially funded by the German Research Foundation (DFG) – 440719683.

OceanRAIN version 1.0, OceanRAIN-W - Water cycle components, 73 along-track parameters for 8 ships, 6.83 million minutes in total, temporally continuous data for each ship, 1-minute-resolution

OceanRAIN version 1.0, OceanRAIN-M - Number Concentration Particle Size Distribution and Precipitation Microphysics, 37 along-track parameters plus 128 size bins for 8 ships, 692.000 precipitation minutes in total, temporally discontinuous data for each ship, 1-minute-resolution

OceanRAIN version 1.0, OceanRAIN-R - ODM470 Raw number count Particle Size Distribution and Precipitation Microphysics, 37 along-track parameters plus 128 size bins for 8 ships, 692.000 precipitation minutes in total, temporally discontinuous data for each ship, 1-minute-resolution